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摘 要 :探讨 构建 足 - 鞋 三 维 复合 有 限 元 模型 的 方法 ,基于 所 建 模型 开展 生物 力学 应 用 研究 ,为 跑鞋 
设计 优化 及 足 部 潜在 损伤 预防 提供 科学 依据 。 基 于 足 - 鞋 三 维 CT 断层 影像 数据 进行 逆向 重建 , 结 
合 实验 三 维 空间 参数 定义 模型 加 载 及 边界 条 件 , 建 立 个 体 化 足 - 鞋 非 完 全 接触 有 限 元 模型 。 细 划 足 
底 和 鞋底 压力 分 区 ,与 实验 测试 结果 对 比 全 方位 验证 模型 的 可 靠 性 ,并 采用 Bland-Altman 分 析 、 
Pearson 相关 性 分 析 以 及 配对 了 检验 进一步 验证 方法 的 一 致 性 及 有 效 性 。 应 用 模型 模拟 分 析 裸 足 / 
着 鞋 状态 ,不同 鞋底 材料 刚度 对 人 体 静 态 站 立 工 况 下 距 骨 区 应 力 的 影响 。 结 果 如 下 :1) 除 前 足 外 ， 
其 他 所 有 分 区 模拟 与 实验 结果 相对 一 致 。 进 一 步 统 计 分 析 发 现 ,两 种 方法 有 显著 相关 性 (r = 
0. 986; P «0.001) ,所 有 压力 差 值 点 均 处 于 95% 一 致 性 界限 区 间 内 , 且 数 据 无 显著 性 差异 (P= 
0. 202) ,表明 方法 的 一 致 性 和 有 效 性 均 相 对 较 好 ;2) 静态 站 立 工 况 下 应 力主 要 分 布 在 第 二 、 三 距 骨 
体 , 其 中 第 三 距 骨 所 受 应 力 最 大 。 与 着 鞋 站 立 相 比 , 裸 足 时 距 骨 区 应 力 峰 值 明 显 上 升 ;3) 随 着 鞋底 
材料 刚度 的 增加 ,又 骨 间 应 力 峰值 差异 也 逐渐 增 大 ,应 力 趋向 集中 于 第 二 、 三 距 骨 。 基 于 上 述 研 究 
结果 ,本 研究 得 出 以 下 结论 :1) 前 足 模拟 压力 不 足 或 BXShREAdx, 后 续 研 究 应 考虑 优化 
约束 条 件 ;2) 裸 足 及 鞋底 刚度 增加 会 导致 距 骨 应 力 异 常 变化 ,建议 在 进行 裸 足 跑鞋 或 碳 板 跑鞋 研 
发 设计 时 ,可 考虑 在 前 足 距 骨 区 域 添加 缓冲 结构 ,协助 优化 路 骨 应 力 分 布 ,降低 潜在 损伤 风险 。 
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Abstract:To develop and validate a 3D comprehensive finite element ( FE) coupled model of the foot- 
shoe complex for revealing their biomechanical interaction and providing scientific basis for running shoe 
design optimization and related foot injury prevention. The collected CT images were used for establishment 
of the subject-specific coupled FE model. Boundary and loading conditions were quantitatively defined 
through 3D gait measurement. Both plantar and shoe sole regions were further divided for model 
verification. The Bland-Altman plot, Pearson correlation analysis, and T test were used for the reliability 
and effectiveness analyses between methods. Based on the established model, the effects of barefoot/shoed 
conditions and shoe sole stiffness on metatarsal stress were further studied. Results are as follows:1) The 
simulated peak pressure showed great consistency with the experimental data except for forefoot regions. 
Statistical analysis found that there was a significant correlation between the two methods (r 20.986; P < 
0.001) ,and that all the pressure difference points were scattered in the 95% limits of agreement , and the 
data had no significant difference ( P 20. 202) , indicating a great agreement between the approaches. 2 ) 
The stress was mainly distributed at the 2™ and 3" metatarsal bones, with the 3" one being the most 
stressed , and it increased in most metatarsal bones while barefoot standing compared to shoed condition. 3) 
As the sole stiffness increased, the peak stress difference between the metatarsal bones increased , and the 


rd 


stress tended to concentrate on the 2" and 3" metatarsal bones. Conclusions are as follows; 1 ) The 
insufficient simulated forefoot pressure may be related to its excessive longitudinal strain, and subsequent 
research should consider further optimization of the constraints. 2) Barefoot and increased sole stiffness 
can lead to abnormal changes in metatarsal stress;thus in terms of running shoe design ,it is suggested to 
consider adding a buffer structure at the metatarsal region for optimizing the metatarsal stress distribution 
and further reducing any potential injury risk. 


Key words :finite element method ; foot-shoe coupled model ; biomechanics ; running shoe ; metatarsal stress 
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足 部 作为 人 体 下 肢 与 外 界 环境 接触 的 始 端 ,在 


洞悉 鞋 具 参数 变化 对 足 - 鞋 内 部 力学 状态 改变 的 影 


日 常 运动 中 发 挥 着 重要 作用 。 同 时 , 足 部 也 是 运 
动 过 程 中 最 常见 的 损伤 部 位 之 一 "” V SERES JE 
部 的 延伸 ,是 运动 过 程 中 足 部 最 直接 、 最 有 效 的 防护 
装备 。 随 着 材料 和 结构 工程 学 的 不 断 发 展 进 步 , 运 
动 鞋 也 从 单一 的 足 部 防护 工具 逐步 成 为 具备 缓冲 减 
震 、 能 量 回 弹 以 及 后 跟 控 制 等 功能 的 综合 运动 装 
备 ”。 一 双 功 能 优异 的 运动 鞋 不 仅 能 够 起 到 足 部 防 
护 作 用 ,还 能 够 显著 地 提升 运动 表现 。 例 如 ,前 期 研 
究 表明 ,Nike Vaporfly 496 跑鞋 在 提升 马拉松 跑 者 跑 
步 经 济 性 上 的 优异 表现 可 能 与 其 在 材料 和 结构 上 的 
创新 有 一 定 关联 。 关 于 运动 鞋 性 能 的 优化 升级 
一 直 以 来 都 是 相关 研究 的 重要 主题 ,其 中 “ 足 - 鞋 - 
地 ”系统 中 的 力学 传递 机 制 以 及 鞋底 材料 选取 与 绥 
冲 回 弹 功效 等 问题 是 主要 研究 热点 |。 

鞋 具 功能 的 每 一 次 优化 升级 都 与 生物 力学 人 研究 
息息相关 ,任何 结构 和 技术 上 的 突破 创新 也 必然 要 
遵循 最 基本 的 生物 力学 机 理 ”] 。 传 统 的 生物 力学 研 
究 方法 虽然 能 够 较为 全 面 地 探究 鞋 具 相 关 特 性 对 足 
部 乃至 下 上 肢 运 动 表现 的 影响 ,但 却 无 法 进一步 
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啊 , 大 大 地 限制 了 足 部 运动 损伤 内 在 机 理 的 挖 气 。 
随 着 有 限 元 法 在 量化 足 部 损伤 疾病 内 在 力学 特征 中 
的 应 用 S ,近年 来 众多 专家 学 者 也 针对 足 - 鞋 交互 
作用 的 有 限 元 模拟 分 析 展 开 了 一 定 探索 ” 。 然 
而 ,目前 足 - 鞋 模型 的 搭建 主要 是 在 足 部 模型 的 基础 
通过 三 维 建 模 所 得 ,并 不 能 较为 真实 地 还 原 运 动 鞋 
的 功能 结构 。 值 得 注意 的 是 ,部 分 研究 在 模拟 过 程 
中 对 模型 进一步 简化 处 理 , 仅 聚 焦 于 鞋底 的 作用 而 
忽略 了 鞋 面 对 足 部 运动 功能 影响 的 重要 性 。 为 解决 
上 述 问题 , 刘 嫌 嫌 应 用 三 维 扫 描 技 术 对 着 鞋 状态 
的 足 部 进行 了 逆向 重建 ,但 是 其 采用 的 足 - 鞋 绑 定 耦 
合 方式 造成 了 鞋 腑 内 部 不 规则 的 凹凸 结构 进而 导致 
其 模拟 效果 偏差 较 大 。 

基于 前 期 探索 ,本 研究 的 首要 目的 在 于 构建 更 
为 真实 完整 的 个 体 化 足 -运动 鞋 三 维 摩擦 耦合 有 限 
元 模型 。 与 此 同时 ,现存 的 足 - 鞋 模型 多 数 仅 采 用 足 
跟 压力 作为 实测 验证 的 单一 参考 指标 , 却 忽 略 了 足 
底 其 他 区 域 以 及 鞋底 区 域 压力 变化 的 验证 。 本 人 研究 
将 对 足 底 和 鞋底 区 域 展 开 进 一 步 划 分 以 更 加 全 面 地 
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验证 模型 的 有 效 性 。 此 外 ,本 研究 将 应 用 模型 模拟 
分 析 裸 足 / 着 圣 状 态 以 及 鞋底 材料 参数 变化 对 人 体 
静态 站 立 工 况 下 踊 骨 区 域 应 力 变 化 影响 , 则 在 为 后 
期 的 足 - 鞋 生物 力学 人 研究 英 定 基础 ,为 运动 鞋 功能 验 
证 提供 可 视 化 平台 。 


1 材料 与 方法 


1.1 研究 对 象 


本 研究 招募 了 1 名 男性 业余 跑 者 (年 龄 :27 岁 ; 
身高 :175 cm; 体 质量 70 kg) 作为 实验 受 试 者 ,其 跑 
姿 为 后 跟着 地 ,有 5 年 户外 /室内 跑 台 跑步 经 验 ,每 
周 跑 量 约 25 km。 受 试 者 无 任何 下 及 畸形 或 心血 管 
疾病 且 测 试 前 半年 内 无 运动 损伤 。 测 试 前 , 受 试 者 
详细 了 解 了 测试 研究 内 容 并 自愿 签署 了 知情 同意 
书 。 本 次 测试 获得 了 作者 所 在 单位 的 科学 伦理 委员 
会 批准 。 


1.2 有限 元 模型 重建 


本 研究 采用 计算 机 断层 扫描 CT 设备 (Optima 
CT540 ,GE Healthcare , Chicago , USA ) 对 受 试 者 右 足 
着 鞋 状 态 下 进行 横断 面 扫描 ,扫描 前 通过 矫正 文具 
使 足 躁 关节 处 于 正中 位 。 将 获得 的 DICOM 格式 影 
像 数 据 导 入 MIMIC 21. 0 软件 (Materialise , Leuven , 
Belgium) 中 进行 足 部 骨骼 、 软 组 织 以 及 运动 鞋 的 三 
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维 几 何 逆向 建 模 。 在 保证 运动 鞋 整体 轮廓 及 其 各 部 
位 厚度 的 基础 上 ,进一步 手动 删除 足 与 鞋 之 间 的 品 
音像 素 ,实现 足 - 鞋 非 完 全 接触 ,建立 有 内 部 腔 隙 的 
运动 鞋 模型 。 此 外 ,为 提升 模拟 运算 效率 ,本 研究 对 
模型 进行 部 分 简化 :第 二 至 第 五 远 端 趾 骨 与 其 相应 
近 端 趾 骨 融合 ;运动 袜 与 软组织 融合 ;运动 鞋 分 为 鞋 
面 和 鞋底 两 部 分 " "9 。 通 过 Geomagic Wrap 2017 软 
件 (3D Systems, South Carolina, USA ) 对 上 述 逆向 建 
模 获 得 的 面 模型 进行 表面 光滑 处 理 并 导入 Solid- 
works 2020 软件 ( Dassault Systèmes , France) 中 做 进 
一 步 实体 化 处 理 。 

为 实现 足 骨 的 相互 连接 和 相对 运动 ,根据 足 中 
解放 结构 对 足 部 主要 软骨 、 币 带 以 及 足 底 筋 膜 进行 
逆向 建 模 ,其 中 软骨 定义 为 实体 模型 而 韧带 和 足 底 
筋 膜 为 线 体 模 型 。 如 图 1 所 示 , 本 研究 的 足 部 有 限 
元 模型 共 由 20 块 足 骨 、20 块 软骨 、1 块 软组织 .66 
FREAR 5. 条 足 底 筋 膜 组 成 。 随 后 ,将 整体 足 - 鞋 
实体 模型 导入 ANSYS Workbench 2021 软件 
( ANSYS Inc. , Canonsburg, USA) 的 Mechanical Model 
TEETH XA S BE RU REI RS IS HP 3 e R 
体 单元 模拟 。 除 支撑 板 采用 六 面体 单元 外 ,其 余 实 


体 模型 均 采 用 四 面体 单元 进行 网 格 处 理 。 其 中 , 软 
2H 2H 3e I EE VL S PCM Ft 70 AJ 5 mm , £r REOS 
3.5 mm ,软骨 为 2 mm ,最 终 模 型 共计 358 322 个 节点 
以 及 208 225 个 网 格 单元 。 


软骨 
(绿色 实体 ) 


gn 
(黄色 短线 体 ) 


C EER 


图 1 E-a ERA A RIRA 


Fig.1  Three-dimensional finite element model of the foot and sports shoe complex 
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1.3 有限 元 模型 材料 属性 


本 研究 足 - 鞋 模型 的 所 有 组 织 模块 均 设 定 为 单 
一 各 向 同性 的 线 弹 性 材料 ,并 通过 杨 氏 模 量 和 泊 松 
比 来 定义 其 材料 属性 。 其 中 , 根据 皮质 骨 与 松 质 骨 
所 占 体 积 比 定义 足 骨 的 材料 力学 参数 ,软骨 、 韧 带 、 
软组织 .运动 鞋 以 及 支 撑 板 的 材料 参数 也 分 别 取 自 
前 期 有 限 元 研究 中 ,具体 见 表 1。 此 外 ,为 进 
一 步 明 确 本 研究 所 需 鞋 底 材料 参数 并 探究 不 同 鞋 底 
材料 对 踊 骨 区 域 应 力 的 影响 ,其 杨 氏 模 量 参 数 在 现 
有 基础 上 作 +10% , +20% 处 理 , 即 2. 490 (0% ) 、 
2.739( +10% ) , 2.241( —1096) ,2. 988( +20% ) , 
1.992 MPa( -2096) 。 

表 1 足 -运动 鞋 有 限 元 模型 材料 参数 
Tab.1 Material properties assigned to each component in 


the foot-sports shoe finite element model 


部 位 杨 氏 模 量 E/ — Madrid 密度 id 
MPa mm (kg * m^?) 
TE 11.76 0.35 - 9 400 
鞋底 2.49 0.35 = 2 300 
骨骼 7300 0.30 一 1 500 
软骨 1 0.40 - 1050 
Wir 260 0.40 18.4 937 
EUR fs 350 0.40 58.6 937 
软组织 1.15 0.49 - 937 
支撑 板 17 000 0.10 过 5 000 


1.4 有 限 元 模型 加 载 及 边界 条 件 


鉴于 模拟 过 程 中 足 底 压力 的 变化 与 足 部 在 三 维 
平面 内 的 相对 角度 具有 一 定 相 关 性 ,本 研究 采用 
8 摄像 头 的 VICON 三 维 动作 捕捉 系统 (Vicon Met- 
rics Ltd. , Oxford , UK) 记 录 受 试 者 静态 站 立时 足 部 反 
光 Mark. 点 的 位 置 参数 ,并 通过 足 部 刚体 坐标 系 
(XYZ) 在 空间 坐标 系 上 的 投影 向 量 进一步 计算 足 部 
在 空间 矢 状 面 和 冠状 面 的 朝向 角 (a = -4.73°;pB = 
-1.15?) P? 。 实 验 测 试 过 程 中 要 求 受 试 者 静态 站 
Nr. SUI H A AOT 7H Fel S , ARAR, 目 视 前 方 。 
如 图 2(a) 所 示 , 足 跟 骨 下 缘 中 心 点 定义 为 刚体 坐标 
系 原点 , 工 轴 指向 第 一 路 骨 与 第 五 趾 骨 连 线 中 间 点 ， 
Z a3 ELT. X Of" CIR] EY lae ECT. XZ 轴 所 在 
平面 。 该 方法 也 将 用 于 后 期 动态 模拟 过 程 中 足 部 转 
动 载荷 边界 条 件 的 确定 。AMITI 三 维 测 力 台 (AMTI， 
Watertown , Massachusetts, USA) 与 VICON 动作 捕捉 
系统 同步 进行 ,用 于 采集 受 试 者 静态 站 立时 垂直 地 
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面 反作用 力 (343. 00 N) ,并 根据 力 的 等 效 互 换 原 理 
将 其 施加 在 支撑 板 底面 中 心 位 置 ,方向 垂直 向 上 。 
前 期 研究 显示 ,人 体 在 双 足 静态 平衡 站 立时 ,下 上 肝 小 
腿 三 头 肌 的 作用 力 约 为 足 部 载荷 的 一 半 ““”。 
此 ,本 研究 于 跟 骨 节点 处 竖 直 向 上 以 集中 载荷 的 形 
式 施加 171.50 N 的 作用 力 用 于 模拟 三 头 肌 作用 力 。 
模型 边界 条 件 设 置 如 图 2(b) 所 示 , 其 中 软 组 
织 \ 胜 骨 以 及 腓骨 的 上 表面 设 定 为 完全 国定 ; 文 撑 板 
设 定 为 仅 可 上 下 移动 ,其 他 方向 则 被 完全 约束 。 软 
组 织 表面 与 鞋 腔 内 以 及 鞋底 与 支撑 板 上 表面 之 间 均 
定义 为 摩擦 接触 ,摩擦 系数 为 0.6 00000, 
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| S) 
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(a) (b) 
到 2 Jis xj EE PROCESAR PARTE 
Fig.2 The application of boundary and loading conditions for 


the foot-sports shoe finite element model 
1.5 ”有限 元 模型 验证 


Pedar 鞋垫 式 压力 传 感 系统 和 Emed 压力 平板 
( Novel GmbH , Munich , Germany ) 分 别 用 于 采集 受 试 
者 静态 站 立时 足 底 和 鞋底 的 压力 分 布 ,测试 频率 均 
设 定 为 100 Hz 并 经 过 实验 室 空 载 标定 ,将 测量 值 与 
有 限 元 模拟 得 到 的 压力 云图 进行 对 比分 析 , 从 而 验 
证 上 述 足 - 鞋 有 限 元 模型 的 有 效 性 。 其 中 ,鞋垫 根据 
足 部 解剖 结构 被 分 为 前 足 内 侧 区 、 前 足 外 侧 区 .中 足 
区 以 及 足 跟 区 ;鞋底 根据 其 结构 特征 被 分 为 前 内 侧 
区 .前 外 侧 区 后 内 侧 区 以 及 后 外 侧 区 ,以 更 加 全 面 
地 将 有 限 元 模拟 值 与 实测 值 匹配 验证 。 此 外 ,本 研 
究 对 实验 测试 和 模拟 分 析 的 压力 值 进行 数据 一 致 性 
分 析 ( Bland-Altman analysis) ,相关 性 分 析 ( Pearson 
correlation analysis) 以 及 配对 样本 工 检验 (Pair-Sam- 
ples T Test) ,以 进一步 评估 有 限 元 模拟 结果 的 可 靠 
性 和 有 效 性 。 其 中 ,采用 MedCalc19. 7. 2 软件 ( Med- 
Calc, Ostend, Belgium ) 进行 一 致 性 分 析 , 采用 
SPSS17.0 软件 (SPSS, Chicago, IL, United States) 3# 
行 相关 性 及 配对 样本 了 检验 分 析 ; 相 关 性 系数 定义 
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为 Ir| 二 0.35 为 弱 相 关 ,0. 36 & Irl <0. 67 为 中 等 
相关 ,0. 68 & Iri «1.00 为 强 相关 '” ; 显著 性 水 平 设 
定 为 P «0.05, 


2 研究 结果 


2.1 模型 验证 结果 


如 图 3 所 示 , 通 过 对 比 有 限 元 模拟 和 实验 测试 
的 压力 云图 发 现 , 两 种 方法 所 测 得 的 足 底 及 土 底 压 
力 分 布 基本 相同 。 其 中 , 足 底 压力 集中 于 足 跟 区 ,其 
次 为 前 足 内 侧 .前 足 外 侧 以 及 中 足 区 ;鞋底 压力 集中 
于 内 侧 区 域 , 即 前 内 侧 区 和 后 内 侧 区 ,其 次 为 前 外 侧 


区 和 后 外 侧 区 。 
MPa 
0.111 00 Max| 
0.103 00 
0.095 50 
0.087 50 
0.079 60 
10.071 60 
(0.063 60 
0.055 70 
0.047 70 
(0.039 80 
0.031 80 40 
0.023 90 
0.015 90 15 
(0.007 96 


(a) 足 底 压力 


MPa 


0.186 0 Max| 


0 Min 


(b) 鞋底 压力 


到 3 区 域 划 分 以 及 有 限 元 模拟 和 实验 测试 压力 云图 对 比 


Fig.3 The divided regions and comparison between predicted 


pressure distribution and the experimental pressure data 
对 比分 析 各 区 域 压 力 峰 值 发 现 , 足 底 中 足 区 和 
后 跟 区 有 限 元 模拟 结果 与 实验 测试 数据 的 相近 度 较 
好 , 且 压 力 峰值 的 相对 误差 均 小 于 10% ;然而 ,前 足 
内 外 侧 有 限 元 模拟 结果 与 实验 测试 数据 的 相近 度 相 
对 较 弱 ,其 中 前 足 内 侧 区 的 压力 峰值 相对 误差 为 
-29.07% ,前 足 外 侧 区 的 压力 峰值 相对 误差 为 
-32. 2196 ,最 终 足 底 压 力 峰 值 平均 误差 为 17. 8996 
(3€ 2) 。 鞋 底 有 限 元 模拟 结果 与 实验 测试 数据 如 表 
3 所 示 ,综合 相对 误差 和 绝对 误差 结果 发 现 ,鞋底 材 
料 参数 为 2.739 MPa ( + 10% ) 时 模拟 结果 与 实验 测 
试 数据 的 相近 度 最 高 。 
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表 2 有 限 元 模拟 和 实验 测试 足 底 各 分 区 压力 峰值 


Tab.2 Predicted peak pressures and the experimental 


pressure data in each plantar region 


实验 测试 / 有 限 元 模拟 / 


名 称 MPa MPa 相对 误差 /% 
足 底 

前 足 内 侧 区 0.086 0.061 -29.07 

前 足 外 侧 区 0. 065 0.044 -32.21 
PEK 0.047 0.043 -8.51 

后 跟 区 0.113 0.111 -1.77 

平均 误差 /% 了 P| =17.89 


表 3 有 限 元 模拟 和 实验 测试 鞋底 各 分 区 压力 峰值 


Tab.3 Predicted peak pressures and the experimental 


pressure data in each sole region 


实验 测试 / 有 限 元 模拟 / 


si 
名 称 MPa MPa 相对 误差 /% 
鞋底 096 

前 内 侧 区 0.182 0.186 2.20 

前 外 侧 区 0.120 0.112 -6.67 

后 内 侧 区 0.189 0.173 -8.47 

后 外 侧 区 0.154 0. 142 -7.79 
平均 误差 /% Y [5| =6.28 
TENE +10% 

前 内 侧 区 0.182 0.198 8.79 

前 外 侧 区 0. 120 0.117 -2.50 

后 内 侧 区 0.189 0.182 -3.70 

后 外 侧 区 0.154 0.150 -2.60 
平均 误差 /% MEZ | -4.39 
鞋底 -10% 

前 内 侧 区 0.182 0.176 -3.30 

前 外 侧 区 0.120 0. 108 — 10.00 

后 内 侧 区 0.189 0. 165 -12.70 

后 外 侧 区 0. 154 0.136 - 11.69 
平均 误差 /% Y P| =9.42 
鞋底 +20% 

前 内 侧 区 0.182 0.207 13.74 

前 外 侧 区 0. 120 0. 120 0.00 

后 内 侧 区 0.189 0. 190 0. 53 

后 外 侧 区 0.154 0.155 0. 65 
平均 误差 /% Ei P| =3.73 
TEE -20% 

前 内 侧 区 0.182 0.166 -8.79 

前 外 侧 区 0. 120 0. 104 -13.33 

后 内 侧 区 0.189 0.157 -16.93 

后 外 侧 区 0.154 0.131 - 14.94 
平均 误差 /% Y P, | =13.49 

上 力学 学 
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对 实验 测试 和 模拟 分 析 的 足 底 及 鞋底 压力 峰值 
做 进一步 统计 分 析 后 发 现 ,两 种 方法 具有 显著 相关 
性 (r=0.986;P<0.001) ,两 种 方法 的 压力 差 值 平均 
值 (0. 006 Mpa) 接 近 于 0, 所 有 压力 差 值 点 均 处 于 
95% 置信 区 间 (limits of agreement) 范围 内 ,日 数据 无 
显著 性 差异 (P=0.202)。 上 述 结果 表明 本 研究 两 
种 方法 的 一 致 性 以 及 模拟 仿真 的 有 效 性 均 较 好 。 


2.2 路 骨 应 力 结果 


图 4(a) ~(b) 所 示 为 裸 足 和 着 鞋 状态 以 及 不 同 
前 过 观察 看 
三 路 骨 


en A 
， 静 态 站 立 工 况 下 应 力 集中 分 布 在 第 二 、 


着 鞋 +10% 


yr 
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体 ， 其 中 第 三 蹈 骨 所 受 的 应 力 最 大 。 与 裸 足 状 态 下 
相 比 ， 着 鞋 站 立时 蹈 骨 区 应 力 峰 值 明 显 下 降 
(2. 876 MPa) ,然而 随 着 鞋底 材料 刚度 的 增加 ,路 骨 
区 应 力 峰值 也 逐渐 上 升 ,但 均 小 于 裸 足 时 应 力 峰 值 
(5.096 MPa), 

此 外 ,通过 进一步 分 析 各 踊 骨 应 力 峰值 发 现 ( 见 
表 4, 图 4c) ,与 裸 足 状态 下 相 比 ,着 鞋 站 立时 各 踢 骨 
应 力 峰 值 差异 较 小 ,然而 随 着 鞋底 材料 刚度 的 增加 ， 
第 一 、 四 以 及 第 五 踊 骨 应 力 峰值 逐渐 下 降 , 第 二 ,三 


骨 应 力 峰 值 逐 渐 上 升 , 路 骨 间 应 力 峰 值 差 异 增 大 ， 
应 力 趋向 集中 , 但 仍 小 于 裸 足 时 跌 骨 间 应 力 峰 值 
差异 。 


mee》 关 性 snas) 视 忆 
mum) 着 畦 -10%,-20% D 着 鞋 +10%,+20% 


—U0OOOOOoooooooo 


图 4 裸 足 / 着 鞋 状 态 下 有 限 元 模拟 分 析 踢 骨 区 等 效应 力 峰 值 及 其 分 布 变化 


Fig.4 Predicted peak metatarsal von Mises stress and its distribution under barefoot or shoed situations 
表 4 ， 裸 足 / 着 鞋 状态 下 有 限 元 分 析 各 路 骨 等 效应 力 峰 值 


Tab.4 Predicted peak von Mises stress of each metatarsal 


under barefoot or shoed situations 


WE AHE — GE FIE FIE AIE 


MPa +0% -10% +10% -20% +20% 
第 一 路 骨 1.133 1.163 1.163 1.160 1.165 1.159 
第 二 踊 骨 3.720 2.543 2.536 2.543 2.529 2.547 
=E 5.096 2.876 2.869 2.876 2.859 2.879 
第 四 踊 骨 3.080 1.332 1.349 1.317 1.369 1.302 
第 五 路 骨 1.385 0.576 0.598 0.559 0.627 0.541 


3 ”讨论 与 分 析 


足 - 鞋 有 限 元 模拟 分 析 的 价值 在 于 它 能 够 较为 
真实 地 揭示 各 种 工 况 下 足 与 鞋 的 生物 力学 相互 作 
用 ,进而 为 鞋 具 材 料 性 能 测试 整体 结构 优化 以 及 是 
部 损伤 预防 提供 科学 依据 “” 。 本 研究 基于 足 - 茜 
实际 CT 扫描 影像 ,对 足 部 骨 骼 、 组 织 以 及 运动 娃 几 
何 外 形 进行 高 度 还 原 , 通 过 三 维 建 模 软件 对 足 部 主 
要 软骨 韧带 以 及 足 底 筋 膜 进行 优化 模拟 ,结合 运动 
生物 力学 三 维 空间 参数 对 模型 加 载 和 边界 条 件 进行 
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量化 定义 ,实现 了 个 体 化 足 -运动 鞋 非 完 全 接触 三 维 
复合 有 限 元 模型 的 构建 。 

模型 验证 是 模型 建立 后 最 重要 的 步 又 之 一 , 通 
过 将 模拟 计算 值 与 实验 测试 值 进 行 对 比 ,分 析 两 种 
方法 的 一 致 性 ,以 期 实现 模拟 仿真 最 大 程度 地 接近 
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果 证 实 了 前 脚掌 着 地 方式 将 显著 增加 趾 骨 区 应 力 。 
李 镜 东 等 "针对 裸 足 状 态 下 足 前 掌 在 不 同 着 地 角 
度 下 跌 骨 区 应 力 变化 的 有 限 元 模拟 分 析 进 一 步 表 
明 , 随 着 前 掌 着 地 角度 增 大 其 中 骨 关 节 应 力 集中 现 
象 将 越发 明显 , 距 骨 应 力 骨折 风险 也 将 越 大 。 本 研 


现实 运动 特征 。 足 跟 区 压力 峰值 是 前 期 足 部 有 限 元 
模型 验证 最 常用 的 参数 指标 '""*  。 本 研究 模拟 足 
跟 区 压力 峰值 与 实验 测量 值 基本 一 致 ,相对 误差 仅 
为 -1.77% ,初步 表明 该 模型 有 效 性 较 好 。 然 而 , 通 
过 对 足 底 和 鞋底 压力 区 域 做 进一步 划分 验证 发 现 ， 
尽管 鞋底 各 区 域 以 及 足 底 中 足 和 后 足 区 相对 误差 均 
较 小 ,前 足 内 外 侧 模 拟 和 测试 压力 值 的 差异 却 相 对 
较 大 ,通过 分 析 推 测 其 主要 原因 如 下 。 其 一 ,本 研究 
采用 ANSYS 静 力 结构 进行 有 限 元 模拟 分 析 , DSL 
腓骨 及 软组织 上 表面 ,通过 在 支撑 板 底面 施加 反 向 
作用 载荷 挤 压 足 - 鞋 模型 进而 模拟 静态 站 立 工 况 。 
然而 ,模拟 过 程 中 前 足 部 位 出 现 的 纵向 应 变 可 能 间 
接 抵消 了 该 区 域 的 部 分 载荷 ,进而 降低 了 模拟 的 足 
底 压力 。 后 期 研究 在 通过 静 力 结构 模拟 足 部 静态 站 
立 工 况 时 应 考虑 在 前 足 区 上 表面 添加 约束 条 件 以 降 
低压 力 误差 。 其 二 ,本 研究 旨 在 建立 非 完全 接触 足 - 
生 有 限 元 模型 ,因此 建 模 过 程 手动 删除 的 足 与 鞋 腔 
之 间 的 噪声 像素 过 多 可 能 会 增 大 内 部 腔 隙 ,导致 内 
部 接触 存在 不 足 并 降低 了 模拟 的 足 底 压力 。 最 后 ， 
基于 前 期 相关 研究 结论 ,本 研究 在 模拟 时 仅 考 虑 小 
腿 三 头 肌 的 作用 力 ,而 忽略 了 其 他 足 部 肌 群 的 影 
响 53 2 。 然 而 ,可 以 预见 , 跟 骨 结 节 处 集中 载荷 的 
增 大 能 够 帮助 减少 前 足 的 纵向 应 变 , 进而 降低 此 区 
域 的 压力 误差 。 尽 管 如 此 ,本 研究 采用 数据 一 致 性 
分 析 , 相关 性 分 析 以 及 配对 样本 T 检 验 进一步 综合 
验证 了 该 足 - 鞋 模型 有 限 元 模拟 结果 的 可 靠 性 和 有 
效 性 。 

中 上 骨 作 为 中 前 足 的 主要 组 成 部 分 之 一 ,对 于 足 
部 的 应 力 传导 以 及 承重 缓冲 有 着 至 关 重 要 的 作 
用 站。 近年 来 , 裸 足 跑步 运动 不 断 兴起 。 与 着 鞋 跑 
相 比 , 裸 足 跑 最 大 的 特点 在 于 其 “回归 自然 ”的 前 掌 
落地 跑步 姿态 。 跑 步 落地 方式 的 转变 能 够 有 效 地 降 
低 裸 足 跑 过 程 中 地 面 的 冲击 波 进而 减少 了 跑步 对 下 
肢 及 足 部 的 损伤 ”1。 然 而 ,前 脚掌 着 地 的 跑步 
姿态 或 将 显著 加 剧 着 地 过 程 中 前 足 距 骨 关节 的 应 
71 ,潜在 增加 该 区 域 损伤 的 风险 i。MORALES- 
ORCAJO 等 针对 裸 足 状态 下 不 同 着 地 方式 对 足 
部 应 力 分 布 的 影响 展开 了 有 限 元 模拟 分 析 , 研究 结 


究 结 果 表 明 ,即使 在 静态 站 立 工 况 下 , 裸 足 的 外 侧 中 
骨 应 力 (第 三 至 第 五 中 骨 ) 也 高 达 着 鞋 状态 下 的 2 
首 。 换 言 之 ,通过 鞋底 的 缓冲 作用 ,运动 畦 能 够 在 一 
定 程 度 上 有 效 降低 路 骨 应 力 ,潜在 降低 其 损伤 风险 。 

目前 专业 跑鞋 鞋底 多 采用 全 掌 碳 板结 构 增 加 其 
TUS RUE , 旨 在 进一步 提高 其 鞋底 回 弹 效率 ,进而 提 
升 跑 步 经 济 性 。 然 而 ,本 研究 发 现 随 着 鞋底 
材料 刚度 的 增加 ,第 二 三 路 骨 应 力 峰 值 逐 渐 上 升 而 
第 一 、 四 以 及 第 五 中 骨 应 力 峰 值 逐 渐 下 降 , 踊 骨 区 域 
应 力 趋 向 集中 。 本 团队 前 期 针对 猎 科 动物 足 掌 的 组 
震 特 性 也 展开 了 有 限 元 模拟 分 析 , 人 研究 结果 发 现 猫 
的 足 热 能 够 帮助 优化 其 趾 骨 区 的 应 力 分 布 ,避免 过 
多 的 应 力 集中 “”” 。 基 于 此 ,本 研究 认为 在 进行 跑 
秆 设计 研发 时 可 考虑 在 前 足 踊 骨 区 域 添 加 “仿生 足 
垫 " 结 构 进 而 优化 路 骨 应 力 分 布 , 或 可 协助 降低 由 于 
长 时 间 重 复 载荷 刺激 而 引发 的 跌 骨 损伤 的 潜在 可 能 
性 。 然 而 ,需要 考虑 的 一 点 是 ,目前 本 研究 仅 对 静态 
站 立 工 况 下 的 足 - 鞋 相互 作用 进行 了 探究 ,下 一 步 应 
将 细 化 该 足 - 鞋 模型 ,完善 足 部 肌肉 载荷 ,对 不 同 鞋 
底 厚 度 ,不 同 跑 速 以 及 裸 足 / 着 鞋 跑步 条 件 下 前 掌 
地 时 踊 骨 应 力 变 化 进行 有 限 元 模拟 分 析 , 进一步 验 
证 论点 并 为 跑鞋 设计 优化 及 距 骨 损伤 预防 提供 精确 
指导 。 


4 结 论 


本 研究 基于 个 体 化 骨骼 组 织 参数 及 步 态 特征 ， 
成 功 构建 足 -运动 鞋 非 完全 接触 三 维 复 合 有 限 元 模 
型 ,并 以 此 对 静态 站 立 工 况 进 行 模拟 仿真 ,探索 不 同 
条 件 下 跑鞋 对 踊 骨 区 域 应 力 变 化 的 影响 。 基 于 研究 
结果 ,提出 建议 如 下 。 

1) 模型 验证 时 前 足 模拟 压力 不 足 可 能 与 其 纵 
向 应 变 过 大 有 关 , 后 续 人 研究 应 考虑 进一步 优化 约束 
条 件 。 

2) 裸 足 及 鞋底 刚度 增加 会 导致 踊 骨 应 力 异常 变 
化 ,建议 未 来 在 进行 裸 足 跑鞋 及 碳 板 跑鞋 设计 研发 
时 ,可 考虑 在 前 足 踊 骨 区 域 适 当 添 加 缓冲 结构 ,使 得 
跑步 着 地 时 踊 骨 应 力 分 布 更 加 均匀 ,或 可 降低 潜在 
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